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Abstract The magnetic field in an accretion disk is 
estimated assuming that all of the angular momen- 
tum within prescribed accretion disk radii is removed 
by a jet. The magnetic field estimated at the base 
of the jet is extrapolated to the blazar emission re- 
gion using a model for a relativistic axisymmetric jet 
combined with some simplifying assumptions based on 
the relativistic nature of the flow. The extrapolated 
magnetic field is compared with estimates based upon 
the synchrotron and inverse Compton emission from 
three blazars, MKN 501, MKN 421 and PKS 2155- 
304. The magnetic fields evaluated from pure syn- 
chrotron self- Compton models are inconsistent with 
the magnetic fields extrapolated in this way. However, 
in two cases inverse Compton models in which a sub- 
stantial part of the soft photon field is generated locally 
agree well, mainly because these models imply magnetic 
field strengths which are closer to being consistent with 
Poynting flux dominated jets. This comparison is based 
on estimating the mass accretion rate from the jet en- 
ergy flux. Further comparisons along these lines will be 
facilitated by independent estimates of the mass accre- 
tion rate in blazars and by more detailed models for jet 
propagation near the black hole. 
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1 Introduction 

We believe that the accretion onto black holes and 
other compact objects is driven by magnetic fields and 
we can also estimate the required magnetic field (see 
§ ??). In blazars, we have independent estimates of 
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magnetic fields reasonably close to the black hole at 
say, ~ 100 — 1, 000 gravitational radii. The aim of this 
paper is to explore the notion that we can relate the 
magnetic field estimates in these two regions by suit- 
ably extrapolating disk magnetic fields to the blazar 
emission zone. We then examine the conditions un- 
der which the extrapolated and estimated fields agree. 
Some of the physical issues which we examine include 
those related to Poynting flux dominated jets, the ac- 
curacy of specific one-zone inverse Compton models for 
blazar emission and the accuracy of estimates of the 
magnetic field strengths in the accretion disk near the 
black hole. 



2 Magnetic fields and jet power 

2.1 Disk magnetic field 

The very existence of a jet or wind in a black hole - ac- 
cretion disk system points to the existence of a magnetic 
str ess additional to the conventi onal r<f> stress identified 
by Shakura and Sunvaev ( 19731 ). (We adopt r, <f> and z 
as conventional cylindrical polar coordinates with z = 
being the disk mid-plane.) iKuncic and Bicknell (|2004l ) 



showed that the channelling of accretion power into ei- 
ther coronal heating or an outflow depends upon the 
existence of a 4>z stress (i0 Z ) + over the surface of the 
disk. In principle, this stress can be important even 
when it is numerically small compared to the r<p stress 
{tr<j})- (t<pz) + acts over the surface are a of the disk, 
whereas (t r <j>) acts over the disk height h so that if 
(tcj, z ) + > {h/r)t r ^ then the dynamical effect of + 
dominates. 

Let M a be the disk mass accretion rate and v$ « Vk 
be the approximately Keplerian mean azimuthal ve- 



locitjQ, then the full equation for angular momentum 
transport is: 
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In principle the mass accretion rate can vary with ra- 
dius r if there is substantial mass flux deposited into the 
wind. However, if the wind and/or jet from the accre- 
tion disk is Poynting flux dominated then dM a /dr « 0. 
Suppose that we further assume that the {t^ z ) + stress 
dominates and that the magnetic stresses dominate the 
kinetic Reynolds stresses (—(pv^Vj)). We then obtain 
the following estimate for the magnetic term (B'.B' Z ), 
viz. 
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where m is the mass accretion rate in units of the Ed- 
dington value A'/Edd = ^GMm p / cut, M is the mass 
of the black hole and r g = GM/c 2 is the gravitational 
radius. Hence we have the following numerical estimate 
for the strength of the magnetic field: 



-B'B'^I 2 = 3.1 x 10 4 M^m 1 ' 2 ( - 



-5/4 



(5) 



where the mass of the black hole is \0 s M%Mq. Co- 
incidentally, this is the s ame dependence of m a gnetic 
flux density assumed by iBlandford and Pavnel |l982h 
in their paper on centrifugally driven winds from accre- 
tion disks. 

2.2 Jet power 

Let us assume that the jet is Poynting flux dominated 
and that the base of the jet outflow is confined to 
radii r m < r < r out . We can determine the jet en- 
ergy flux Fe by integrating the mean Poynting flux 
(S z ) w —^{B'^B'^Vfj, over the cross-sectional area of 



1 We utilize the mass weighted averaging of the MHD equa- 
tions in which direct means are denoted by a bar and mass- 
weighted means by a tilde. Turbulent components are denoted 
by a prime and ensemble averages by angular brackets. See 
iKuncic and Bicknelll d2004h for details. 



the jet utilising equation |4} for (B^B^). Hence, 
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This equation can thus be used to estimate the mass 
accretion rate given the jet energy flux. Then elimi- 
nating M a from equation |4]) we have, for the magnetic 
field as a function of cylindrical radius ro at the base of 
the jet: 
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The determination of the base magnetic field in 
terms of the jet can be used in conjunction with a blazar 
emission model (which gives the jet energy flux) to de- 
termine the magnetic field in the jet extrapolated from 
the accretion disk. This is a valid procedure since the 
jet energy flux is conserved irrespective of entrainment 
in the region b e tween the disk and the blazar emission 



region lBicknelll (|1994I ). This is important since the par- 



ticles in the initially Poynting dominated jet may be 
entrained en-route. 



3 Relativistic magnetised jet flow 

In order to connect the magnetic field at the base of 
the jet to the field further out in the blazar emission 
region we need to consider the evolution of the field 
along the Poynting flux dominated jet. In order to do 
this we consid er the equatio ns of relativistic axisymmet- 
ric MHD (see Mestell (1999)) with a newtonian gravita- 
tional field. Let fl be the angular velocity of the jet, p 
its rest mass density, v p = (v r , 0, v z ) the poloidal veloc- 
ity, Bp = (B r ,0,B z ) the polodial magnetic field and T 
the Lorentz factor. Then we have the following stream 
line constants: 



pvp 
Bp 

Ratio of poloidal rest mass flux 

to poloidal magnetic flux (8) 
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rB^ 
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Rate of angular momentum transport 
per unit poloidal flux tube (9) 
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Angular velocity of field lines 
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One of the important surfaces in such a flow is the 
Alfven surface where the poloidal gas velocity equals 
the poloidal Alfven speed (v p = va)- The Alfven sur- 
face is also equivalently characterised by AirTrj 2 / p = 1 
and — /3/4-7T = arjTr 2 . Denoting the values of dynami- 
cal variables at the Alfven surface by a subscript A, the 
azimuthal magnetic field (the dominant perpendicular 
component) is given by: 



(1 - V p A lT AP ) 



(11) 



The asymptotic form of this expression is obtained for 
r > r A , p < PA and T ~ T A , giving 



rB^ » -AirarjTApr 2 /pa 



(12) 



The value of rB^ at the disk surface is obtained assum- 
ing that r <C ta and p 3> pa giving 



roB^fi w -4irarjT A r A 



Hence, 
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The other streamline parameter that we require is the 
ratio of azimuthal to vertical field at the disk surface. 
The first step is to evaluate the streamline constant 



a = \l - =\l — 

pr 



(15) 



Assuming that B^ is at most of order B p and that the 
initial poloidal velocity is much less than the Keplerian 
velocity (v p <C i>k,o = r ty we have 



a w CI 
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Then, dividing equation (TT^|) by tqB z ^ and using 
B z , = r]~ x pov Z fi gives 
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In both equations (|T3|) and (|14j) we have the similar 
terms pr / p A T A and p A v ai~a/poVz,oTq which we can 
treat approximately as follows. Let the equation of a 
poloidal stream line by r = R(z). The continuity equa- 
tion tells us that 



d.z 



R(z) 



p(r, z)v z (r, z)r dr = 
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If we assume that over a jet cross-section the density 
and vertical velocity are constant, i.e. p(r, z) w p(z) 
and v z (r, z) w v z (z), then 



Thus along any poloidal streamline, pv z r 2 
Consequently 
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Furthermore, we assume that by the time the flow has 
reached the Alfven surface it is already relativistic al- 
though possibly not with a high Lorentz factor. Hence, 
we assume that va ~ c. Also, in the asymptotic regime 
where T ^> 1, v z rs c to an even better approximation. 
Therefore, in the asymptotic region 
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The ratio of the azimuthal and vertical components of 
magnetic field in the disk is: 



B, 
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Given this estimate of the ratio of toroidal to vertical 
field components at the base of the jet we can use the 
estimate for (B ( f }t0 B Z fi) + on the accretion disk surface 
to determine the asymptotic azimuthal field using equa- 
tion ||2"2"1). The result is: 
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where tq refers to the value of r at the base of the 
streamline. In order to compare this estimate of mag- 
netic field with the one-zone models which invoke a sin- 
gle magnetic field strength we define an area weighted 
flux density 
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where i?j ct is the radius of the jet in the blazar emis- 
sion region. In order to carry out this integral one needs 
to integrate^o/rg) -1 / 2 with respect to the streamline 
radius r across the jet in the blazar region; this in 
turn requires that we know the relationship between 
the streamline radius and the radius at the base of the 
streamline, i.e. r(ro). We make the simplifying as- 
sumption that the expansion is linear. This should not 
be a great source of error since the integrand only de- 
pends weakly on tq. With this assumption and with 



p(z)v z (z)R 2 (z) w constant 



(19) 
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x = rin/fout as the ratio of inner to outer jet radii: 
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The function /(x) varies by a factor of 2 from its mean 
value « 0.24 for 0.019 < x < 0.94 so that the precise 
value of x is unimportant. However, for illustration we 
consider two different values in the following. 

The magnetic field estimate in equation (|26p repre- 
sents the magnetic field in the laboratory frame. For 
comparison with blazar estimates the magnetic field in 
the plasma rest frame is obtained by dividing this esti- 
mate by the Lorentz factor T. 



4 Comparison with models of blazar emission 

Our aim is to compare the estimate of the magnetic 
field using the above theory with the values derived 
from inverse Compton models of blazar high energy 
emission. Inverse Compton emission results from the 
scattering of "soft" photons by high energy electrons 
and one may distinguish three qualitatively different 
models for blazar emission: 

1. Exernal inverse Compton (EIC) models in which 
the soft photons originate from outside the region 
of high energy electrons. 

2. Synchroton Self Compton (SSC) models in which 
the soft photons are produced by the same pop- 
ulation of high energy electrons which scatter 
them. 

3. Local Inverse Compton (LIC) models in which the 
soft photons originate from a region which is adja- 
cent to the region of high energy electrons. How- 
ever, the regions may be physically related. For 
example, one can envisage a region of jet contain- 
ing embedded shocks. The high energy electrons 
would then be in a sub- volume of a larger volume 
which is the source of soft photons. 

In the literature the distnction between LIC emis- 
sion and EIC emission is not usually made. How- 
ever, we think that the distinction is useful in view 
of the close spatial and physical relationship that 
may exist between the source of soft photons and the 
scattering region. Examples of LIC emission mod- 
els include the decelerating flow mo dels proposed by 
Kazanas and Georganopoulosl (120061) and the 'blob in 
jet' model of iKatarzvhski. Sol, and Kud (|200lh . 



Table [T] summarises the comparison of blazar and 
extrapolated accretion disk estimates for three blazars. 
Different groups have independently estimated blazar 
parameters using different models (SSC and LIC) for 
MKN 501. In the comparisons we have used two dif- 
ferent values of x = ri n /r out ; however, as indicated ear- 
lier this does not make great deal of difference to the 
extrapolated magnetic field estimates. We have also 
indicated which of the models are SSC and which are 
LIC; the last column of the table indicates the ratio of 
the extrapolated magnetic field to the magnetic field 
estimated from the particular blazar model. 



5 Discussion 

One of the obvious conclusions from Table Q] is that 
the SSC models do not fare very well in estimating the 
magnetic field in the blazar region. On the other hand 
the two LIC estimates of magnetic flux density agree 
well with the extrapolated magnetic field. Whilst such 
an agreement is welcome it is not especially surprising. 
The estimate of the magnetic field in equation (|27jl is 
based on a Poynting flux dominated jet which remains 
so in the asymptotic regime - co nsistent with the orig- 
inal n on-relativistic solution of iBlandford and Pavnel 
1 19821 ). It is well known that pure SSC models yield low 
magnetic energy densities and high particle energy den- 
sities that are incompatible with an important Poynting 
flux. The LIC models on the other hand effectively re- 
duce the particle energy density and increase the mag- 
netic energy density as a consequence of the larger flux 
of soft photons. The LIC models therefore imply com- 
parable values of magnetic flux density and energy den- 
sity, which are more compatible with a substantial com- 
ponent of Poynting flux. Hence, the agreement with the 
extrapolated field is not surprising. In order to provide 
a more stringent test of the extrapolation technique de- 
veloped here it is desirable to have an independent es- 
timate of the accretion rate (rather than one based on 
the energy flux) which can be used to estimate the base 
magnetic field (using equation (jl])) independently of the 
energy flux. This could be provided, for example, by 
the X-ray and/or optical luminosity associated with the 
accretion disk. 

On the other hand, the reasonable comparison be- 
tween extrapolated and estimated magnetic fields bodes 
well for a continuation of this approach in which sev- 
eral of the simplifying assumptions can be relaxed. In 
particular the gravitational field of the black hole can 
be incorporated into a fully general relativistic treat- 
ment, a small mass flux can be incorporated and the 
consistent development of the model for the flow near 
the black hole may also be included. 
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Source 


Blazar 


-Rjct 


Tjct 


K 


^IC 


F E 


^in/^*out 


-Bext 


Bext/Bic 




Model 


10 16 cm 




cm -3 


C 


10 43 ergs s _1 




c 




MKN 501 


ssc 1 


0.53 


10 


4.3 x 10 3 


0.06 


7.0 


0.1 


1.3 


22 


MKN 501 


ssc 1 


0.53 


10 


4.3 x 10 3 


0.06 


7.0 


0.5 


1.7 


28 


MKN 501 


LIC 2 


1.0 


7 


100 


0.15 


0.33 


0.1 


0.19 


1.3 


MKN 501 


LIC 2 


1.0 


7 


100 


0.15 


0.33 


0.5 


0.25 


1.7 


MKN 421 


SSC 3 


1.0 


10 


4 x 10 4 


0.06 


100 


0.1 


2.4 


39 


MKN 421 


SSC 3 


1.0 


10 


4 x 10 4 


0.06 


100 


0.5 


3.0 


50 


PKS 2155-304 


LIC 4 


0.65 


25 


160 


0.12 


0.68 


0.1 


0.12 


0.8 


PKS 2155-304 


Lie 4 


0.65 


25 


160 


0.15 


0.68 


0.5 


0.15 


1.0 



Table 1 Comparison of the extrapolated magnetic field i? cx t and the inverse-Compton estimated magnetic field 
Bic for three blazars. The jet radius, Lorentz factor, number density and B\q are provided by the cited models; 
the jet energy flux is inferred from those models. The extrapolated magnetic field is calculated using the above 
theory and the rati o of extrapolated to inverse Compto n magnetic field is given in the last col umn of the ta- 



ble. R efer ences: 1: Bicknell. Wagner, and Groves! ( 2001 ); 2: Katarzvhski. Sol, and Kus ( 2001 ); 3 Maraschi et al 



(jl999h : 4: Uharonian et all |2005l) Model 1 



